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ABSTRACT 


Finding the first generation of stars formed out of pristine gas in the early Universe, known as Population III (PopIII) stars, is one of the most 
important goals of modern astrophysics. Recent models suggest that PopIII stars may form in pockets of pristine gas in the halo of more evolved 
CY) galaxies. Here we present NIRSpec-IFU and NIRSpec-MSA observations of the region around GN-z11, an exceptionally luminous galaxy at 
z = 10.6, which reveal a >So detection of a feature consistent with being HeIIA1640 emission at the redshift of GN-z11. The very high equivalent 
M width of the putative Hell emission in this clump (170 A), and the lack of metal lines, can be explained in terms of photoionisation by PoplII 
<L stars, while photoionisation by Popll stars is inconsistent with the data. It would also indicate that the putative PopIII stars likely have a top-heavy 
initial mass function (IMF), with an upper cutoff reaching at least 500 Mo. The PopIII bolometric luminosity inferred from the Hell line would be 

e ~ 2x 107? Lo, which (with a top-heavy IMF) would imply a total stellar mass formed in the burst of ~ 6 x 10° Mo. We find that photoionisation by 
the Active Galactic Nucleus (AGN) in GN-z11 cannot account for the Hell luminosity observed in the clump, but can potentially be responsible 

for additional Hell emission observed closer to GN-z11. We also consider the possibility of in-situ photoionisation by an accreting Direct Collapse 

I Black Hole (DCBH) hosted by the Hell clump; we find that this scenario is less favoured, but it remains a possible alternative interpretation. We 
© also report the detection of a Lya halo stemming out of GN-z11 and extending out to ~2 kpc, as well as resolved, funnel-shaped CII] emission, 


> likely tracing the ionisation cone of the AGN. 


Jun 2023 
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> 
C^) 1. Introduction lent widths (EW(HeIL41640) > 20 A), unaccompanied by metal 
un . i lines (Tumlinson & Shull 2000, Oh et al. 2001, Tumlinson et al. 
ON The formation of the first stars and galaxies marks a fundamental 2001, Schaerer 2003, Nakajima & Maiolino 2022). However, 
© transitional phase in cosmic history, during which the Universe other models have also predicted that primordial gas might frag- 
© evolved from a relatively simple state into the highly structured ment more efficiently, resulting into the formation of cooler and 


XO System we observe today (Bromm & Larson 2004). Almost all less massive stars (e.g. Clark et al. 2011; Greif et al. 2011, 2012; 
© baryonic matter formed after the Big Bang and the “dark ages" is Stacy et al. 2016) 


CNA composed of hydrogen and helium. Consequently, the very first 
CN stars to form must have condensed out of these elements, pri- 
* marily as a result of cooling due to the small amounts of H, that 
«= could form at such early epochs (e.g. Kashlinsky & Rees 1983). 
< These stars, called Population III (PopIID), are potentially very 
$} massive, up to several 100 Mo, as a consequence of inefficient : > : 
C3 gas cooling and hence poor fragmentation of pre-stellar cores at have eluded detection so far, with some claims not confirmed 
these early epochs (e.g. Bromm et al. 1999, Abel et al. 2002,Tan by subsequent studies (Sobral et al. 2015, Bowler et al. 2017, 
& McKee 2004; Hirano et al. 2014; Susa et al. 2014). Indirect Harikane et al. 2018). However, the unprecedented sensitivity 
studies of Pop III stars via stellar archaeology confirm these pre- of JWST is pushing the frontier of observations to high redshift 
dictions, suggesting that their masses extend up to 1000 Mo with and very faint galaxies, discovering very low metallicity systems 
a characteristic mass > 1 Mọ (Rossi et al. 2021; Pagnini et al. (Vanzella et al. 2023), indicating that the detection of Poplll stars 
2023). Such massive PopllI stars are expected to be very hot May be within reach of this observatory (see also e.g. Zackrisson 

and leading to an energetic photoionising spectrum, capable of ®t al. 2011). 
doubly ionising helium. Indeed, an expected signature of PoplIII Further encouraging expectations about the detectability of 
stars is prominent Hell nebular lines, with very large equiva- —PopllI stars have been obtained by recent cosmological simula- 
tions. Indeed, two main formation phases of PopllI stars have 
rm665 @cam.ac.uk been proposed: firstly, PoplII stars could have formed in dark 


While observations have found potential chemical finger- 
prints of the enrichment produced by the first generation of stars 
(Beers & Christlieb 2005; Frebel et al. 2007; Cooke et al. 2011; 
Hartwig et al. 2015, 2018; Salvadori et al. 2023; Saccardi et al. 
2023; Klessen & Glover 2023), PoplIII stars in early galaxies 


* 
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matter mini haloes at redshifts z~10-30 (e.g. Tegmark et al. 
1997; Bromm et al. 2009; Pallottini et al. 2014; Jeon et al. 2015, 
Vikaeus et al. 2022; Ventura et al. 2023; Trinca et al. 2023b). 
Secondly, PoplII stars could still have formed at later times, 
down to redshifts of about z ~ 3-6, in pristine/low-metallicity 
pockets of gas due to the highly inhomogeneous nature of metal 
enrichment in the Universe, as well as a consequence of re- 
accretion of pristine gas (e.g. Tornatore et al. 2007a; Johnson 
2010; Liu & Bromm 2020, Sarmento et al. 2018; Jaacks et al. 
2019, Katz et al. 2022, Volonteri et al. 2023). In these cases, 
PoplII stars are expected to be found in the haloes of more mas- 
sive galaxies, and their hard ionizing spectrum may be ioniz- 
ing mildly enriched gas, resulting in weak emission of metal 
lines (Nakajima & Maiolino 2022, Katz et al. 2022, Jaacks et al. 
2018). 


Although in the latter scenario PopllI stars are expected to 
be present even at relatively low redshifts (possibly down to 
z~3 ; Cassata et al. 2013, Liu & Bromm 2020), the probabil- 
ity for a halo to host PoplII stars is higher at very high redshift 
(Trinca et al. 2023a). Within this context, GN-z11 is an excep- 
tional galaxy at z=10.6 which offers an optimal environment 
for searching PoplIII stars. Identified initially by Bouwens et al. 
(2010); Oesch et al. (2016), GN-z11 is the most luminous galaxy 
at z>10 in all HST fields (including all CANDELS and Fron- 
tier Fields, Oesch et al. 2018). Recent JWST-NIRCam imaging 
from the JWST Advanced Deep Extragalactic Survey (JADES; 
D. Eisenstein, in prep.) revealed that GN-z11 is characterized 
by an unresolved central source, a more extended component 
(200 pc in radius) with Sérsic index n=1, and a more diffuse 
component extending for ~0.5’ to the NE (the ‘haze’), which 
may not be associated with GN-z11 and could simply be a fore- 
ground system (Tacchella et al. 2023). Bunker et al. (2023) ob- 
tained low and medium resolution spectroscopy of GN-z11 with 
JWST-NIRSpec in the multi-object spectroscopy (MOS) mode 
using the micro-shutter assembly (MSA), revealing a spectrum 
with a wealth of nebular emission lines in various ionization 
stages, and also characterized by some unusual line ratios. In- 
terestingly, the spectrum reveals Lya emission, extended along 
the slit for about 0.2" (~ 0.8 kpc) towards the SW (the posi- 
tion angle of the shutter was 19.9°). Additional NIRSpec MSA 
observations targeting GN-z11 have led to the identification of 
an active galactic nucleus (AGN) in GN-z11 (Maiolino et al. 
2023), primarily through the detection of high ionization lines 
(NeIV), semi-forbidden lines tracing very high densities consis- 
tent with the Broad Line Regions (BLR) of AGN, and the detec- 
tion of the CII* 41335 fluorescent emission, typically seen only 
in AGN. The spectrum of GN-z11 has also revealed the presence 
of a powerful outflow, traced by the redshifted (~ 530 kms !) 
resonantly scattered emission of both Lya and CIVA1550, and 
blueshifted absorption of CIV (~ 800 — 1000 km s*!). The sec- 
ond NIRSpec MSA observation had a position angle nearly per- 
pendicular to the first observation, and did not show obvious ev- 
idence for Lya extension on either side of GN-z11. 


A more detailed analysis of the 2D spectrum of GN-z11 in 
the first observation revealed the tentative emission of a line lo- 
cated at ~0.5” to the NE of GN-z11, at a wavelength consis- 
tent with HeIL11640 at the same redshift of GN-z11 (see Sect. 
3.1). This intriguing finding prompted a DDT proposal to obtain 
NIRSpec-IFU spectroscopy of GN-z11 with both the G140M 
(3.3 h) and G235M (10.6 h) gratings. The IFU datasets deliver 
superior spatial information relative to the MSA (although the 
sensitivity is lower than the MSA, because the IFU has about 
50% lower throughput). 
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In this paper we report a preliminary analysis of the new 
NIRSpec IFU data of GN-z11, which confirm the HelL41640 
detection at z=10.6, as well as additional interesting features as- 
sociated with the Lya and CIII] morphology. 


Foreground 
galaxy 
(z=2.028) 


MSA O 


pointing 
accuracy 


Fig. 1: F277W NIRCam image of the central ~ 1.7" around 
GN-z11 illustrating the location of the MSA shutters during the 
first NIRSpec-MOS observations in February 2023 (tilted boxes, 
with boxes of the same color being associated with the same con- 
figuration), extraction region of the Hell clump (red square) and 
of the Hell extended emission (yellow box). The white circle in- 
dicates the absolute uncertainty of the MSA target acquisition 
(hence the shutters’ footprints are uncertain by this amount). 
North is up and East is to the left. The vertical dashed orange 
line marks the boundary beyond which fewer dithers from the 
NIRSpec IFU observations are available because of the guid- 
ing problem reported in Sect. 2. Note that the full IFU field of 
view is larger than shown here (>3”x 3” taking into account the 
dithers). 


Throughout this work, we use the AB magnitude system and 
assume a Planck flat ACDM cosmology with Qm = 0.315 and 
Ho = 67.4 km/s/Mpc (Planck Collaboration et al. 2020). With 
this cosmology 1" corresponds to a transverse distance of 4.08 
proper kpc at z — 10.6. 


2. Observations, data processing and data analysis 


GN-z11 was observed with NIRSpec IFU (Jakobsen et al. 2022; 
Boker et al. 2022) on May 22"4 — 23" 2023 under the DDT 
program 4426 (PI: Roberto Maiolino). The observational setup 
utilised a medium cycling pattern of 10 dithers for a total in- 
tegration time on source of 10.6 h with the medium-resolution 
grating/filter pair G235M/F170LP, and 3.3 h with the medium- 
resolution grating/filter pair G140M/F100LP. These configura- 
tions cover the spectral ranges of 0.97—1.89 and 1.66-3.17 um, 
with a nominal spectral resolution of R ~ 1000. The detector 
was set with the improved reference sampling and subtraction 
pattern (IRS*), which reduces significantly the read out noise 
with respect to the conventional method (Rauscher et al. 2012). 
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Fig. 2: Top panel: Signal-to-noise 2D prism spectrum of GN-z11 from the February 2023 MSA observation, reprocessed to use 
only the exposures in which GN-z11 is in the upper and lower shutters (hence discarding the exposure with GN-z11 in the central 
shutter), so to avoid self-subtraction of extended signal. Bottom panel: 1D spectrum extracted from the region delimited by the two 
dashed horizontal lines, at about 0.45" from GN-z11 in the NE direction, illustrating the tentative detection (3.7 o integrated) of a 
feature at 1.903 um. The red vertical line indicates the location of HelI41640 expected at the redshift of GN-z11, z=10.603 (which 
is undistinguishable from the redshift of 10.600 inferred for the grating spectrum). 


We also checked that no bright foreground sources were at the 
location of the MSA quadrants where they could significantly 
contaminate the IFS spectra. 

Although JWST has an excellent pointing accuracy (i.e. 
0.1", 1 c radial! Rigby et al. 2023), and the telescope was com- 
manded to the location of GN-z11 provided by the NIRCam and 
HST imaging and GAIA-referenced astrometry, the observations 
were affected by a guiding problem of the observatory, because 
the guide star turned out to be a binary (report from the FGS 
team); this introduced a large off-set of about 1.4" towards the 
SE. This offset located GN-z11 close to the edge of the IFU FoV. 
This restricted the zone around GN-z11 that could be explored, 
lowering the S/N of regions which were covered by fewer dither 
positions, or excluding completely some regions from the FoV. 
More specifically, out of 10 dither positions, GN-z11 itself was 
located right at the edge of the FoV in three of them, and in an- 
other four of them GN-z11 is so close to the edge that the SW 
extension of the Lya was mostly outside the FoV. Although in 
this paper we attempt an analysis of the SW extension based on 
the incomplete data, this will need to be validated by planned 
future observations that will repeat the dithers that were more 
severely affected (within the context of the same program 4426). 
As a consequence, the main focus of this paper is on the extended 
emission towards the NE. 

Because of the guiding problem, the world coordinates given 
by the pipeline in the header of the final IFU cube are off by 
about 1.4". Therefore, we corrected the world coordinates by re- 
aligning the cube to the NIRCam image, by collapsing the cube 


! https://jwst-docs.stsci.edu/ 
jwst-observatory-characteristics/ 
jwst-pointing-performance 


in the same band as in the NIRCam filter (although the centroid 
of the continuum does not change with wavelength). 

Raw data files were downloaded from the MAST archive 
and subsequently processed with the JWST Science Calibration 
pipeline? version 1.8.2 under CRDS context jwst_1068.pmap. 
We made several modifications to the default reduction steps 
to increase data quality, which are described in detail by Perna 
et al. (2023) and which we briefly summarize here. The default 
pipeline is occasionally subject to over-subtraction of elongated 
cosmic ray artefacts, through the default circular correction ap- 
plied to *snowball-like' cosmic ray hits. To avoid this problem, 
we patched the pipeline to fit ellipses to all flagged regions con- 
sisting of five or more adjacent pixels; regions with best-fit el- 
lipses having axial ratios smaller than 0.1 are removed from 
the list of *snowball-like' cosmic rays. Count-rate frames were 
corrected for 1/f noise through a polynomial fit. Outliers were 
flagged on the individual 2D exposures, using an edge-detection 
algorithm similar to LAcosMic (van Dokkum 2001). We calculate 
the derivative of the count-rate maps along the x-axis of the de- 
tector (on the scale of a few pixels, the x-axis is a good approx- 
imation for the dispersion direction). The derivative was nor- 
malised by the local flux (or by 3 times the rms noise, whichever 
was highest), and we rejected the 95'^ percentile of the result- 
ing distribution — the percentiles are calculated over the entire 
detector area (see F. D' Eugenio et al. 2023 for details). In addi- 
tion, we made the following corrections to the *cal.fits files after 
Stage 2: we mask pixels at the edge of the slices (two pixels 
wide) to conservatively exclude pixels with unreliable sflat cor- 
rections. We also mask regions that are affected through leak- 


? https://jwst-pipeline.readthedocs.io/en/stable/ 
jwst/introduction.html 
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age from failed open MSA shutters. Note that the data quality 
flag '*MSA FAILED OPEN ' fails to account for all failed open 
MSA shutters. At the same time, several regions marked as be- 
ing affected by failed open shutters do not show any impact from 
leakage, and we include those regions to increase S/N in the fi- 
nal combined cube. Finally, we mask regions that are affected 
by persistence in all dither frames, e.g. from cosmic rays. The 
final cubes were combined using the ‘drizzle’ method with pixel 
scales of 0.06" and 0.10", for which we used an official patch to 
correct a known bug?. 

We also re-analyse the NIRSpec MOS observations of GN- 
z11 obtained in February 2023 within the JADES survey, but 
with focus on the 2D extraction along the shutters. These obser- 
vations were described in Bunker et al. (2023). Here we only re- 
mind that those observations were obtained with three configura- 
tions of the MSA (Jakobsen et al. 2022; Ferruit et al. 2022; Bóker 
et al. 2023), each with 3-shutters nodding. The low-resolution 
prism and the three medium-resolution gratings were used, for 
a total observing time of 6.9 hours with the prism and 3.45 
hours with each of the gratings. In the following, to explore 
the 2D information we use only the prism, as the gratings ex- 
posures are shallower. The processing of the MSA data is also 
described in Bunker et al. (2023), here we only recall that we 
have used the pipeline developed by the ESA NIRSpec Science 
Operations Team and the NIRSpec GTO Team. We also recall 
that the spatial pixel size in the MSA 2D spectra is 0.1" and 
that each shutter is 0.2" wide. The standard background sub- 
traction method does not allow to explore extended emission 
along the shutters. Indeed, the nodding technique moves the tar- 
get by one shutter position, i.e. by 0."53 along the shutters di- 
rection. When reciprocally subtracting two exposures in which 
the source is located in adjacent shutters, this results into a deep 
negative trace of the subtracted source spectrum at +0.5” from 
the (positive) spectrum of the source, wide at least a much as 
the PSF and its wings (~0.2”) or more, depending on the source 
size. This makes the standard 2D spectrum unusable to explore 
any extended emission beyond ~0.2” from the source. There- 
fore, in order to explore the 2D information along the shutter, 
and avoid self-subtraction of the extended emission, in this pa- 
per we use a different background subtraction process relative 
to standard process adepted in Bunker et al. (2023). Specifically, 
we do not use the exposure with the source in the central shutter, 
and use only the exposures in which the source is in the upper 
and lower shutters (i.e. the most distant ones) to perform mutual 
background subtraction. This ensures that the subtracted, nega- 
tive traces of the GN-z11 spectrum, are located at ~ 1" from 
the (positive) spectrum, and enables the study of spectral fea- 
tures along the shutter up to about 0.6" from the source. This 
method obviously slightly reduces the sensitivity of the data, as 
we use only 2/3 of the data and also the noise in the background 
is increased by a factor of V2, but is far superior relative to the 
standard method in exploring extended emission. 

The tilted boxes in Fig.1 show the footprint of the shutters in 
the four MSA configurations of the MOS observation in Febru- 
ary 2023, overlaid on the F277W image of GN-z11 (Tacchella 
et al. 2023). The green circle indicates the target acquisition un- 
certainty of the MSA, meaning that the footprints indicated in 
the figure may actually be displaced in any direction by that 
amount. The dashed vertical orange line indicates the region be- 
yond which fewer exposures from the NIRSpec IFU DDT pro- 
gramme are available (and at the edge of the field of view) be- 
cause of the telescope guiding problem discussed above. 


5 https: //github.com/spacetelescope/jwst/pull/7306 
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When extracting spectra and emission line maps, we re- 
move background and continuum emission with the following 
procedure. We extract the background from an area of the IFU 
free from artefacts and by sigma-clipping outliers. The resulting 
spectrum is smoothed with a rolling median of 25 spectral pix- 
els. This background is normalized to the continuum in the vicin- 
ity of the spectral feature and subtracted. Since the background 
has some residual structure across the IFU field of view, also 
because of the possible continuum contribution from individual 
sources (especially in the case of GN-z11), we further optimize 
the continuum subtraction by linearly interpolating the contin- 
uum within a range of a few 0.01 um around the line of interest. 
In the case of the Lya the latter step may be more problematic on 
GN-z11, as the galaxy continuum is asymmetric relative to Lya 
due to IGM absorption (which also affects the continuum red- 
wards of the Lya because of the damping wing); however, we are 
primarily interested in the extended emission, where the galaxy 
continuum is not detected. We have further explored the subtrac- 
tion of the continuum around Lya by fitting the continuum only 
on the red side of the line, and the result does not change, except 
for a slight increase in the noise on the extended component, as 
expected from the poorer sampling of the background. 


3. Results 
3.1. Hell detection in the halo of GN-z11 


We start by illustrating the tentative detection of HelIA1640 in 
the 2D spectrum along the shutters in the NIRSpec MSA obser- 
vation of February 2023 where we show the S/N 2D spectrum 
of GN-z11 in the top panel of Fig. 2. As already discussed, this 
was reprocessed to use only the exposures in which GN-zl11 is 
in the upper and lower shutters (hence discarding the exposure 
with GN-z11 in the central shutter), to prevent self-subtraction 
of extended emission. The 2D spectrum shows a tentative sig- 
nal at 1.905 um, which consistent with the expected location 
of HeIIA1640 at the redshift of GN-z11 (z=10.603), peaking at 
about 0.5" from GN-z11. The bottom panel shows the spectrum 
extracted from the region of the 2D spectrum delimited by the 
two horizontal dashed lines in the upper panel. The black dot- 
ted line indicates the noise level per spectral pixel. The extracted 
spectrum confirms the emission feature at 1.905 um. The inte- 
grated significance of the feature is 3.70, as reported in Table 2. 
We note that this feature is totally impossible to detect with the 
standard background subtraction procedure adopted in Bunker 
et al. (2023), as in that case the deep negative trace of the GN- 
z11 subtracted spectrum falls exactly at this location: this is the 
first negative trace, next to the (positive) spectrum of GN-z11, in 
the top panel of Fig.1 of Bunker et al. (2023); with our method 
we avoid the negative traces nearest to the GN-z11 spectrum (i.e. 
those at +0.5” from the spectrum of GN-z11 in the 2D spectrum 
of Bunker et al. (2023)). We also note that the absolute flux and 
spatial centroid of the feature are affected by some additional un- 
certainty associated with the 2D extraction. Specifically, the fea- 
ture is at the edge of the region beyond which self-subtraction 
can still happen even by adopting our 2-shutters methodology; 
this can result into both flux uncertainty and the centroid being 
slightly artificially shifted towards GN-z11. Secondly, in some 
exposures the feature falls next or beneath the shutter bar (see 
Fig. 1), this can also affect the flux and centroid. Finally, the 
aperture losses applied to the spectrum are based on the point- 
like assumption and by using the location of GN-z11; hence if 
the Hell emission is slightly extended this may result in addi- 
tional flux uncertainties. 
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Fig. 3: Top panels: Maps of the HelIL41640, Lya, and CIII]41900 lines (see text for details). The red star marks the position of the 
continuum of GN-z11. Contours are at 2, 3 and 4 ø for the Hell and Lya maps, and at 2, 4, 6 and 8 o for the CHI] map (the rms is 
estimated from the statistics in the IFU field of view, outside this region). Note that most of these spaxels are independent, hence the 
significance of the integrated emission of each feature is significantly higher. The two white boxes indicate the extraction apertures 
of the two spectra shown in Fig. 4 and marked in red and yellow in Fig. 1. The vertical dashed segment indicates the region beyond 
which there are fewer exposures because of the guiding problem discussed in Sect. 2. Bottom panels: overlay of the contours from 
the emission line maps onto the F277W NIRCam map from Tacchella et al. (2023). 


Fig. 3 top-left shows the maps at the expected wave- 
lengths of Hell, Lya and CII] obtained from the NIRSpec-IFU 
observation. Specifically, Hell was extracted from 1.9003um 
to 1.9035um, CIM] from 2.2099um to 2.2173um, Lya from 
1.4115um to 1.4152um (the extraction of CII] and Lya is on a 
slightly broader range as the former is a doublet, hence slightly 
broader although still unresolved at our resolution, and Lya has 
a broad profile). Note that the Lya map has been centered on 
the Lya offset observed in GN-z11 (-550 km s^!). The maps are 
also slightly smoothed with a Gaussian Kernel with a width of 
1.5 pixels (0.09") . 


Note that, since these features are extended, the significance 
of the integrated emission of the individual clumps is higher than 
that given by the contours. 


The Hell map shows a plume extending for about 0.25" to- 
wards the West of GN-z11 (red star). This is the area with fewer 
exposures (due to the guiding problem), so it should be treated 
with great care. Yet, if confirmed, such Hell plume would extend 
well beyond the size of the galaxy identified by Tacchella et al. 
(2023) and may be tracing gas photoionized by the AGN. The 
weak CIII] (and no CIV emission seen either at this location), 
would indicate that this is very low metallicity gas photoionized 
by the AGN. 


The most intriguing feature is the clump located at about 
0.6" (2.4 kpc) to the NE. This is very close to the location of 
the emission serendipitously found along the shutter of the MSA 
prism observation (Fig. 2), especially taking also into account 
the MSA target acquisition uncertainty and the fact that the MSA 
observation still suffers signal self-subtraction at locations far- 
ther than 0.5" from GN-z11 (as well the additional uncertanties 
discussed in the previous section). We have extracted the spec- 
trum from this region, by taking a 0.24" x 0.24" square aper- 
ture (red box in Fig. 1, ‘Hell clump extraction’, smaller white 
box in Fig. 3 and Tab. 1). The resulting spectrum is shown in 
the top-left panel of Fig. 4, showing the detection of a line at 
Aops = 1.902 um, whose integrated significance is 5.5 c (Table 
2). This corresponds to the wavelength of HelIA1640 at a red- 
shift of 10.600, fully consistent with the redshift of GN-z11. 

It is unlikely that this feature is Hæ or [OIII]45007 of a lower 
redshift interloper, specifically at z=1.90 or z=2.80, as other 
bright lines expected at these redshifts are not seen. The ‘haze’ 
identified in imaging by Tacchella et al. (2023) is probably a 
lower redshift interloper, but is offset (closer to GN-z11) rela- 
tive to the putative Hell clump; we indeed identify an emission 
line at 2.330 um peaking on the ‘haze’ that we tentatively iden- 
tify as [OIII] at z23.65 (see Appendix B for details), consistent 
with the photometric redshift inferred by Tacchella et al. (2023). 
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Fig. 4: Spectra extracted from the apertures indicated in Fig. 1, zoomed around the expected wavelength of Hell at the redshift of 
GN-z11. Specifically, the left panel shows the spectrum extracted from the aperture centered on the “Hell clump’ (red box in Fig. 1 
and smaller white box in Fig.3); the right panel shows the spectrum extracted from the larger aperture encompassing the extended 
Hell emission in the NE quadrant of GN-z11 (yellow box in Fig.l and larger white box in Fig.3). The vertical red bars indicate 
the expected wavelength of Hell at the redshift of GN-z11, z=10.603, and the redshift corresponding to the centroid of the line in 
the Hell clump, z=10.600. The dotted lines indicate the noise per spectral pixel in each spectrum inferred from the rms directly 


measured in the spectrum around the Hell line. 


The large foreground galaxy at z=2.28 to the NW of the clump 
does not have any line emission at this wavelength (see spec- 
trum in Appendix B), and any weak emission feature potentially 
associated with that galaxy would also show much stronger Ha 
associated with it, as well as many other stronger emission lines. 

We have also explored whether some Lya emission is leak- 
ing out of the IGM absoprtion at this location. We found a ten- 
tative detection of Lya, but which requires additional data for 
confirmation, so we do not discuss this any further. 

It is difficult to compare the fluxes of the putative Hell seen in 
IFU clump and the detection in the MSA. In terms of wavelength 
they are fully consistent. In terms of flux it is more difficult. At 
face value, the fluxes reported in Tab.2 give a flux in the MSA 
lower (with large error), as expected by the fact that the MSA 
shutter only picks part of the Hell clump. However, because of 
the uncertainty on the exact location of the MSA shutter (see 
Sect.2 and white circle in Fig.1) it is not known what is the ac- 
tual fraction of the putative Hell clump that entered the shutters. 
Moreover, the potential absolute flux calibration issues plaguing 
the MSA, as discussed above, introduces additional uncertain- 
ties that prevent a proper comparison. Finally, we note that the 
MSA-Prism observation is more sensitive than the IFU observa- 
tion (as discussed in Sect. 1) hence it may be detecting additional 
Hell emission outside the clump, not seen in the IFU cube. 

In addition to the clump at ~0.6” to the NE from GN-z11, the 
Hell emission appears to have additional, more extended emis- 
sion, as seen in the map (Fig. 3 top-left), mostly distributed the 
NE quadrant of GN-z11, possibly with a fainter, less significant 
clump, about 0.3" south of the Hell primary clump and about 
0.4” east of GN-z11. We have extracted the spectrum from this 
larger area by using the large aperture marked with the yellow 
box in Fig. 1 (larger white box in Fig.3). The resulting spec- 
trum at the wavelength of Hell at z=10.600 is shown in the right 
panel of Fig. 4. In this case the putative Hell emission is detected 
at 6.10 and with a total flux that is about two times higher than 
in the clump. We note that the putative Hell emission detected 
in this larger aperture is certainly not polluted by Hell emission 
in GN-z11. To begin with, the Hell emission of GN-z11 is very 
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faint (Maiolino et al. 2023; Bunker et al. 2023); contamination 
by GN-z11 would imply that all other strong lines seen in GN- 
z11 would be detected, and much brighter than Hell. Secondly, 
in Appendix C we also extract a spectrum by modifying this 
large aperture so to remove the region closer than 0.3" (more 
than six times the PSF radius at this wavelength) from GN-z11, 
and the detection remains unchanged (Fig.C.1), therefore con- 
firming that the weak Hell emission from GN-z11 does not con- 
tribute at all to the spectrum extracted in this region. 


3.2. Lya halo and morphology 


Although not the focus of this paper, in middle panel of Fig. 3 
we also show the map of the Lyo emission. Here we simply col- 
lapse the background/continuum-subtracted channels around the 
redshifted peak observed in GN-z11 (3.1), while a more in depth 
analysis of the full Lya profile will be presented in a forthcoming 
paper. The Lya emission is clearly extended, partly in the West- 
ern direction relative to GN-z11 (red star), but mostly towards 
the SW, as expected from the MOS data when taken along this 
orientation (Bunker et al. 2023). The extension towards the SW 
is very sharp and could potentially trace an accreting filament. 
However, as discussed in the next section, the fact that we also 
see lop-sided CIII] emission in this direction favors the interpre- 
tation that the emission in the SW stems from gas illuminated in 
the ionization cone of the AGN in GN-z11. We however warn 
that this region falls in the area with less exposure because of 
the telescope guiding problem discussed in Sect. 2, so the emis- 
sion in this region could potentially extend farther to the SW. A 
proper investigation requires additional data that will be obtained 
in the repeated observations. 


At a low surface brightness level Lya extends smoothly and 
over a larger area also towards the North, with a peak at about 
0.4" to the NE of GN-z11. The nature of the extended Lya emis- 
sion will be discussed more extensively in a dedicated paper. 
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3.3. CIII] extension and the ionization cone 


The CIII]1909 doublet is one of the brightest metal lines ob- 
served in the spectrum of GN-zl11, and is the only one that 
is clearly resolved in the IFU observations. Its continuum- 
subtracted map is shown in the right panel of Fig. 3. The CIII] 
emission is clearly elongated in the NE-SW direction, and most 
prominently extended towards the SW with a funnel-shaped ge- 
ometry (note that the elongation is not an artefact of the data 
processing, as demonstrated in the Appendix where we com- 
pare it with the continuum map extracted around the same wave- 
length). The most plausible interpretation is that this emission 
is tracing the ionisation cone of the AGN in GN-z11. Although 
no CIV is clearly detected in this region (which however might 
be self-absorbed; Maiolino et al. 2023), the emission line prop- 
erties are fully consistent with an AGN photoionisation scenario 
with low ionisation parameter, likley due to the distance from 
the AGN (Feltre et al. 2016; Nakajima & Maiolino 2022). Once 
again, since this extended emission falls in the region where the 
exposure time is shorter because of the telescope guiding prob- 
lem discussed above, therefore the investigation of this extended 
emission will require the additional data that will be obtained in 
the repeated observations. 
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Fig. 5: Observed HeIIA1640 luminosity of the “Hell clump’ and 
its projected distance from GN-z11, compared to the expected 
Hell luminosity as a function of distance from GN-z11, in case 
of AGN photoionization. In this calculation, we use the upper 
limit for the size of the Hell clump and assume a broad range 
of AGN ionising shapes powering the photo-ionisation. Further 
conservative assumptions are that the ‘Hell clump’ is not made 
of unresolved clumps and that the projected distance is the actual 
separation. Clearly, the Hell luminosity observed in the ‘Hell 
clump’ is completely inconsistent with being photoionized by 
the AGN in GN-z11. 


4. The case for pristine gas and Population Ill 


As discussed in Sect. 3.1, the most likely identification of the 
emission line observed at 1.902 um in the halo of GN-z11 is 
HellA1640 emission at z=10.600. In the vicinity of GN-z11, 
specifically the Western plume and some of the extended emis- 
sion towards the NE quadrant, the Hel emission may likely trace 
gas photoionized by the AGN in GN-z11. The fact that no or lit- 
tle emission from metal lines (in particular CIV, CIII] and OIII]) 


is detected in these regions suggests that such circumgalactic gas 
has low metallicity. However, it is also possible that the low in- 
tensity of all metal lines is a consequence of low ionisation pa- 
rameter (Feltre et al. 2016; Nakajima & Maiolino 2022). Yet, 
empirically, when extended Hell emission is observed in other 
haloes surrounding AGN at lower redshift, this is always accom- 
panied by even brighter CIV emission and possibly CII] emis- 
sion (Guo et al. 2020; Fossati et al. 2021). Therefore, the finding 
of extended Hell emission around an AGN at z>10 not accom- 
panied by CIV emission suggests that this is primarily due to 
much lower gas metallicity in the halo. 

However, although the gas closer to GN-z11 may be pho- 
toionized by the AGN, the Hell emission seen in the NE clump 
at ~0.6” (2.4 kpc) from GN-z11 cannot be explained in terms of 
photoionisation by the AGN, simply based on ionising photons 
budget arguments. Indeed, HeII41640 is a simple recombination 
line, which therefore essentially is a counter of photons more en- 
ergetic than 54 eV of the impinging radiation. As a consequence, 
knowing the luminosity of the AGN (Maiolino et al. 2023), as 
well as the distance and size of the cloud (or its upper limit), 
one can estimates the number of Hel ionising photons reaching 
the cloud and the resulting Hell luminosity. We have assumed 
a broad range of AGN ionising spectra (Nakajima & Maiolino 
2022) and that the observed, projected distance gives the actual 
distance of the clump to GN-z11. This obviously is a conserva- 
tive assumption as the actual distance is larger and resulting into 
a lower flux of ionising photons. We also assume that, within the 
clump, the covering factor of the absorbing gas is unity, i.e. that 
the gas in the cloud is not ‘clumpy’, hence all radiation reach- 
ing it is fully absorbed. The result of this calculation is shown in 
Fig. 5 where the blue shaded area shows the expected Hell lu- 
minosity as a function of distance from GN-z11, assuming pho- 
toionisation by the AGN in GN-z11 and by adopting all AGN 
ionising continuum shapes considered in Nakajima & Maiolino 
(2022). We have also assumed isotropic emission, which is an- 
other conservative assumption given that the AGN in GN-z11 is 
type 1 (Maiolino et al. 2023). The luminosity observed in the 
Hell clump (L(HelI) = 7.7 x 10^! erg/s, based on the flux re- 
ported in Table 2) is shown with the golden symbol, which is 
clearly inconsistent with AGN photoionisation by about two or- 
ders of magnitude. Note that we have assumed stable emission, 
but it is very unlikely that the luminosity of the AGN can change 
by two orders of magnitude, even because the variability ob- 
served so far is constrained to the level of 10% (Maiolino et al. 
2023). 

The Hell emission in the clump must therefore result from 
in-situ photoionisation. In Fig. 6 we explore the source of ionis- 
ing photons and gas metallicity by using diagnostic diagrams 
in which the equivalent width of the Hell line is compared 
with various line ratios, specifically: HeII/Hy, CIVA1550/Hell, 
CIII]41909/HeII, OIII|41665/Hell. We compared these quanti- 
ties with the photoionisation models presented in Nakajima & 
Maiolino (2022). Specifically, the purple stars illustrate the case 
of photoionisation by PoplI stars (through metal poor models 
from the BPASS libraries), with a broad range of ionisation pa- 
rameters and densities and with (absolute) metallicities ranging 
from 1075 to 107°. In these models the gas and stars are assumed 
to have the same metallicity. The blue symbols show the mod- 
els for PopIII ionising continua. Dark blue circles are the case 
of PopIII embedded in pristine gas. Light blue squares and cyan 
triangles show the case of PopIII photoionising slightly enriched 
gas, as indicated in the legend. In the case of PopIII, Nakajima & 
Maiolino (2022) use a Salpeter IMF with three different stellar 
ranges: 1-100 Mo, 1-500 Mo, and 50-500 Mg; these are coded 
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Fig. 6: Diagnostic diagrams showing the equivalent width of HeIIA1640 versus various line ratios: Hell/Hy, CIVA1550/Hell, 


CIII]1909/HeII, OIII|41665/Hell. Purple stars show PopII models 


with different (absolute) metallicities (color coded with different 


purple shadings), and where the gas and stars are assumed to have the same metallicity. Blue circles are models of PopIII stars in 
a pristine gas; light blue squares are models of PoplII stars in a mildly enriched gas, specifically Z;4, = 1075; light blue triangles 
are models of PoplII stars in a more enriched gas, specifically Zgas = 107^. The PoplIII markers have different sizes depending on 
the stellar mass range assumed for the IMF, as indicated in the legend. Orange pentagons show the predictions for accreting Direct 
Collapse Black Holes (DCBH) embedded in gas with different (absolute) metallicities. The observational constraints on the Hell 


clump in the halo of GN-z11 are shown with a golden circle. 


with PopIII symbols of different sizes, as indicated in the legend. 
We note that, despite the high mass range explored in Nakajima 
& Maiolino (2022), a Salpeter IMF may not be appropriate for 
PoplII (e.g. de Bennassuti et al. 2014, 2017; Stacy et al. 2016; 
Venditti et al. 2023), and the EW(Hell) for PopIII may therefore 
have values even higher than those obtained by those models. 


The constraints on the Hell clump in the halo of GN-z11 are 
shown with a golden circle. The continuum for estimating the 
EW(Hell) is inferred by extracting and interpolating the NIR- 
Cam photometries in the F150W, F200W and F277W bands (Ta- 
ble 3). Very little flux is observed in the aperture used to extract 
the Hell clump emission, and it is not even clear if such weak 
emission is associated with the Hell clump or with some emis- 
sion from the foreground galaxy to the NE of the clump, or to the 
*haze'. Therefore, the inferred EW is more likely a lower limit. 
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For the Hell/Hy ratio we use the data from the prism, as Hy is 
not covered by the new IFU observations (note that the Hell/Hy 
is a lower limit given that Hy is not detected, see Table 2); in 
this case we measure the EW separately by extracting the pho- 
tometry at the corresponding aperture along the MSA shutter, as 
reported in Table 1. The resulting equivalent widths are given in 
Table 2. The EW(Hell) inferred for the clump and for the MSA 
spectrum are consistent with each other given the uncertainties. 
As said, these EW are likely lower limits. 


Fig. 6 illustrates that the high EW(Hell41640) = 172 + 76 A 
observed in the Hell clump is inconsistent, by orders of mag- 
nitude, with the PopII models, while consistent with the PopIII 
models. Furthermore, it requires a very top-heavy IMF, with an 
upper mass cutoff of at least 500 Mo. Possibly an even higher 
stellar mass cutoff is needed, or, in addition to the high mass cut- 
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off, an IMF more top-heavy than Salpeter, or perhaps a log-flat 
IMF (Chon et al. 2021), are required to explain the very large 
EW(Hell). 

The lower limit on the Hell/Hy ratio (Fig. 6 top left; see 
also Table 2) is nicely consistent with the constraint from the 
EW (Hell), i.e. requiring photoionisation from PopllI stars, and 
possibly with a top-heavy IMF. In contrast, the upper limits from 
the metal lines (see Table 2) are not very constraining in terms of 
metallicity and we cannot exclude that the PoplII stars are em- 
bedded is some moderately enriched gas, although the CIV/Hell 
upper limit is at the borderline with the models of PopllI living 
in a gas with metallicity of 107. 

By using the PopIII models with the two top-heavy IMFs (1- 
500 Mo and 50-500 Mo) we infer the ratio between Hell lumi- 
nosity and total, bolometric luminosity (~ 107? for the two most 
top-heavy IMF models which are better matching the observa- 
tions). This gives a bolometric luminosity of 7.7 x 10? erg s^! 
for the PopIII in the clump. Assuming an Eddington limited lu- 
minosity, this gives a mass of only 6 x 10° Mo. Alternatively, 
by using the ZAMS mass to luminosity conversions for PopIII 
given by Schaerer (2002), the inferred bolometric luminosity 
would give a total stellar mass of 7 x 10° Mo, similar to the 
value obtained above assuming an Eddington limited luminos- 
ity. These values for the total PopIII stellar mass are not far from 
the PopllI stellar mass expected from some recent simulations. 
Indeed, these expect a halo hosting a Popll stellar mass simi- 
lar to GN-z11 (^ 8 x 108 Mo in the extended component, Tac- 
chella et al. 2023), to harbour a fractional mass of PoplII stars 
of about 3 x 107^ — 1 x 1073, i.e. a total PoplII stellar mass of 
~ 2-8x10? Mo (Yajima et al. 2022; Venditti et al. 2023; A. Ven- 
ditti, priv. comm. ). 

It is worth noting that, in addition to the extended Hell emis- 
sion, although the map is noisy, the Hell clump itself seems re- 
solved and possibly extended over ~0.25”, i.e. ~1 kpc. This is 
much larger than the expected size of a single PopIII sub-halo 
and it is much more likely that the extension is associated with a 
cluster of PopHI sub-halos. Indeed, according to recent simula- 
tions, PoplII sub-halos are expected to cluster in regions of a few 
kpc (Tornatore et al. 2007b; Venditti et al. 2023). Alternatively, 
the extension may trace a large clump of pristine gas which is 
being photoionised by a single halo of PopIII stars embedded in 
1t. 

We finally also consider the possibility that the HelI clump 
is ionized in-situ by an accreting Direct Collapse Black Hole 
(DCBH) seed in a pristine, or very low metallicity gas cloud. 
Nakajima & Maiolino (2022) have explored this case in their 
photoionization calculations, and these models are shown with 
orange pentagon symbols in Fig. 6. The darker orange sym- 
bols show the case of DCBH photoionizing pristine (zero metal- 
licity) gas, while lighter shades of orange indicate the case of 
moderately enriched gas, as indicated in the legend. As for the 
PoplII and PoplII case, models explore a broad range of densities 
and ionization parameters. Additionally, models explore a broad 
range of plausible ionizing shapes for the accreting DCBH (we 
defer to Nakajima & Maiolino 2022, for a more detailed scr- 
ription of the models). The DCBH models reach Hell equivalent 
widths much larger than PoplI stars, however not as high as the 
PoplII with the most top-heavy IMF. In particular, the DCBH 
models do not reach the EW (Hell) required to explain the obser- 
vation of the Hell clump. However, the difference is not huge (a 
factor of ~3—5 within the errorbars), and we do not exclude that 
more extreme DCBH models could be consistent with the obser- 
vations. Therefore, although the DCBH scenario is less favoured, 
it remains a possible alternative interpretation. 


5. Conclusions 


We have presented new NIRSpec-IFU observations of GN-z11,a 
remarkably luminous galaxy at z=10.6. We have also presented a 
re-analysis of the 2D MSA spectra of GN-z11 from the JADES 
survey, by employing a data processing technique that enables 
the analysis of extended emission along the shutters. Unfortu- 
nately, a telescope guiding problem during the IFU observations 
resulted in GN-z11 being offset to the western edge of the FoV 
in various dither positions, so the focus of the paper has been on 
the analysis of the extension in the eastern part, while the proper- 
ties of the western extension were only briefly analysed (a more 
detailed analysis of the western extension is deferred to after the 
repeated observations of the affected dithers, within the context 
of the same observing programme). 
The main observational findings are the following: 


— We detect a spectral feature at 1.902 um in a clump at ~ 
0.5 — 0.6" (~ 2 kpc) to the NE of GN-z11. This wavelength 
corresponds to HeIL41640 at z=10.600, fully consistent with 
the redshift of GN-z11. The line is detected both in the MSA- 
prism (3.70) and in the IFU-grating (5.50). 

— The Hell emission is also detected over a more extended 
area in the NE quadrant of GN-z11, possibly with a sec- 
ond (fainter) clump. There is also a Hell plume extending 
towards the West of GN-z11. 

— Lya and CIII] are also clearly extended and well resolved. 
Lya has a sharp feature extending towards the SW of GN- 
z11, where bright CIIT] also extends in a funnel-shaped ge- 
ometry (however, we remind that these extended features to- 
wards the West will need to be re-examined with new data). 
At lower surface brightness Lya also extends over a larger 
region of up to ~2 kpc to the NE of GN-z11, and with a peak 
at about 1.5 kpc from the galaxy. 


We investigate in detail these features, mostly focusing on 
the Hell emission, and reach the following conclusions: 


— While Hell emission close to GN-z11 may be associated 
with photoionization by the AGN, simple photon budget ar- 
guments exclude that this is the case for the Hell clump, leav- 
ing in-situ ionisation as the only possibility. 

— We show that the very high EW(Hell41640) = 172 + 76 A 
and large Hell/Hy ratio cannot be explained in terms of pho- 
toionization by PopII (metal-poor) stars, while they are con- 
sistent with photoionization by PopllI stars. 

— The very high EW(Hell) suggests that the putative PoplIII 
stars must have a top-heavy IMF reaching an upper mass 
cutoff of at least 500 Mo. The non-detection of Hy (hence 
the lower limit on HelI/Hy) also supports this scenario. 

— We infer that the mass of PoplII stars formed in the burst is 
about ~ 6 x 10? Mo, which is not far from the fraction of 
PopIII stars in massive halos at these redshifts expected by 
some simulations. 

— We also consider the alternative possibility of photoioniza- 
tion by a Direct Collapse Black Hole in the Hell clump. This 
scenario is less favoured, as it predicts a lower EW, but not 
by a large factor. Hence this scenario remains a possible in- 
terpretation. 

— Finally, we suggest that the funnel-shaped CIII] extension, 
accompanied by the Lya elongated feature towards the SW, 
traces the ionization cone of the AGN hosted in GN-z11. 
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Table 1: Extraction apertures (given the IFU observation guiding problem, the IFU mask coordinates apply only once the cubes are 
registered to have the GN-z11 continuum matching the absolute coordinates given in Tacchella et al. (2023) 


RA DEC size(’’) PA 
Hell clump 12:36:25.514 +62:14:31.80 0.24x0.24 0 
Hell large ap. 12:36:25.514 +62:14:31.62 0.48x0.60 0 
MSA offset ap. 12:36:25.471 +62:14:31.78 0.20x0.20 19.9 


Table 2: Measured fluxes, equivalent widths and upper limits in the apertures considered in this paper 


Emission line Flux EW 
107?? erg s^! cm? A 
Hell clump small aperture 
Hell41640 5.0+0.94 172+76 
CIV 1550 «34 
OIII]1665 <21 
CIII] 1909 <3.0 
Hell large aperture 

Hell41640 14.1+2.3 63233 
CIV1550 « 8.7 
OIII]1665 « 6.9 
CIII]1909 «1/1 

MSA offset aperture 
Hell41640 3.4+0.9 98452 
Hy <1.37 


Table 3: Continuum fluxes in Hell clump from NIRCam 


Filter Flux [nJy] 
Hell clump small aperture 
F150W 4.4+1.8 
F200W 2.2+1.6 
F277W 3.1+40.9 
Hell large aperture 
F150W 19.7+6.2 
F200W 25.9+8.1 
F277W 25.9+8.1 
MSA offset aperture 
F150W 2.2+1.8 
F200W 3.6+1.6 
F277W 4.9+0.9 
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Appendix A: NIRSpec-ICU contiuum map 


Fig.A.1 shows the map of the continuum obtained by collapsing 
continuum channels near the CIII]1909 emission, with contours 
of the continuum-subtracted CII] line map overlaid. The contin- 
uum map does not show any elongation, demonstrating that the 
CIII] elongation is not an artefact of the data processing. 
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Fig. A.1: Map of the GN-z11 continuum around the CIII] emis- 
sion line with the CIII] contours overlaid (red). Contours are 2, 
4, 6 and 8 c levels. The CHI] emission shows an extension to the 
East while the continuum appears like a point source. 


Appendix B: Interlopers and other emission-line 
sources 


We searched the region around GN-z11 for possible interlop- 
ers. There is a disc galaxy with bright Ha emission 1.1 arcsec 
North-East of GN-z11. We integrate the flux of this galaxy (dark 
red contours in Fig. B.1, left) and show the resulting spectrum 
(top right panel); at the redshift of z — 2.028, there is no strong 
emission line that may contaminate the Hell emission (i.e., near 
A = 1.90 um). Besides, at the location of the Hell ‘blob’ (blue 
contours), the strength of any emission line due to the interloper 
would be even weaker than what shown in the extracted spec- 
trum. Viceversa, if the putative Hell emission in the clump were 
due to some faint line associated with the galaxy at z=2.028, then 
the spectrum extracted from the same aperture of the Hell clump 
would also show much stronger emission of Ha and the several 
other strong emission lines, which are not seen. We can therefore 
exclude that the line we identify as Hell is due to misidentified 
faint line emission at z — 2.028. 

The red contours show the extraction box around the *haze' 
(cf. Tacchella et al. 2023). The corresponding spectrum shows 
an emission line at 2 = 2.331 um. We tentatively identify this 
line as either [OIII]45007 — typically the strongest emission 
line in low-mass galaxies. This would mean the ‘haze’ was a 
z = 3.665 interloper — consistent with the photometric redshift 
of Tacchella et al. (2023). If the line was Ha, we would expect 
to see [OIII]45007 at bluer wavelength (low-mass galaxies have 
typically little dust), but no suitable emission is found. Either 


Article number, page 12 of 13 


way, there is no strong emission line from this interloper near 
the wavelength of Hell. 

Finally, we detect an emission line to the SW of GN-z11 
(bottom panel of Fig. B.1). This source has no continuum and 
therefore no photometric redshift. We are thus unable to identify 
its redshift. Clearly, however, no contamination can arise from 
this region over to the Hell blob. 


Appendix C: Spectrum extracted from large 
aperture excluding regions closer to GN-z11 


The Hell emission in GN-z11 is so faint that it cannot contribute 
to the spectrum shown in the right panel of Fig.4 extracted from 
the large aperture dranw with a yellow box in Fig.1 (Tab.1), or 
else many other much stronger lines from GN-z11 would be seen 
in the spectrum. However, in order to further exclude any con- 
tamination, we have re-extracted the spectrum from the same 
aperture by excluding the region closer than 0.3" (i.e. 6 times 
the radius of the PSF at this wavelength) from GN-z11. The re- 
sulting spectrum is shown in Fig.C.1, which is essentially un- 
changed with respect to Fig.4-right; if anything, the Hell detec- 
tion is slightly more significant in this spectrum. Therefore, it is 
confirmed that the Hell flux from GN-z11 does not contribute to 
the flux observed in that aperture. 
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Fig. B.1: Map of the continuum in a region centred at the ‘haze’, highlighting a bright disc galaxy at z = 2, the haze itself, and a 
region South-West of GN-z11. The spectra of these three sources/regions are shown in the right panels, highlighting the location of 
our Hell detection, as well as other detected emission lines. The galaxy at z — 2 does not have any bright feature at wavelengths 
close to 1.902 um. The haze presents a single emission line which we tentatively identify as [OIII]A5007 at z = 3.66 — which 
would be consistent with the photometric redshift from Tacchella et al. (2023). Finally, we detect an emission line in the SW lobe 
of GN-z11. At the redshift of GN-z11, this line would correspond to a rest-frame wavelength of 0.1684 um — which we do not trace 


to any notable emission line. 
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Fig. C.1: Spectrum extracted from the same extracted aperture at 
the spectrum in the right panel of Fig.4, i.e. the aperture marked 
in yellow in Fig.l (Tab.1), but where the regions closer than 
0.3" to GN-z11 were excluded. The spectrum is essentially un- 
changed with respect to the spectrum extracted from the entire 
aperture. 
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